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\~^ . Abstract 

> 

^\ , A lattice Boltzmann scheme simulating the dynamics of shell 

f^ ' models of turbulence is developed. The influence of high order kinetic 

^Z^ ' modes (ghosts) on the dissipative properties of turbulence dynamics 

>0 , is studied. It is analytically found that when ghost fields relax on 

^^~^ • the same time scale as the hydrodynamic ones, their major effect is 

a net enhancement of the fluid viscosity. The bare fluid viscosity is 
recovered by letting ghost fields evolve on a much longer time scale. 
Analytical results are borne out by high-resolution numerical simu- 
lations. These simulations indicate that the hydrodynamic manifold 
(Ih ' is very robust towards large fiuctuations of non hydrodynamic fields. 

s^ ' Key words: Turbulence, Kmetic Theory, Lattice Boltzmann, Shell 

• ^H , Models. 

X 

H , 

5^ , 1 Introduction 

In the recent years it has been speculated that modern developments in 
(discrete) kinetic theory might yield a new angle of attack to the problem 
of turbulence modeling 0|2]. The rationale behind this idea is as follows. 
All turbulence models based on the Navier-Stokes hydrodynamic equations 
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work on the assumption of a scale separation between the resolved and 
unresolved eddies j^lSJEj. By resolved eddies, we mean excitations at a 
scale larger than the grid size of the simulation. Such an assumption is 
never satisfied in turbulent flows, particularly close to solid boundaries, 
where turbulence production and removal are strongly unbalanced. This 
is the reason why all eddy-viscosity models fail to reproduce accurately 
turbulent statistics in strongly non-homogeneous situations. Kinetic theory 
is at a vantage point to describe such strongly out of equilibrium conditions, 
since it does not require any scale separation between the 'fast' and 'slow' 
degrees of freedom. We wish to emphasize that we are referring to an 
effective kinetic theory dealing with the dynamics of quasiparticles, i.e. 
effective degrees of freedom of the turbulent flow ^^ |^ |22] . 

In any eddy- viscosity model the effect of small eddies on the large ones 
results in a typical diffusion process, only with a much enhanced turbulent 
diffusivity. In contrast to eddy-viscosity, the kinetic approach is centered 
upon the more general and fundamental notion of relaxation, controlled 
usually by a single characteristic time, r entering in the Boltzmann equa- 
tion. The ratio between the relaxation time and the typical hydrodynam- 
ical time, Tfi, is called the Knudsen number, Kn = r/rh- In the limit of 
small Kn, the Boltzmann equation converges to the Navier-Stokes equa- 
tions I^Oll^SI- The major question is: how to derive a suitable effective 
kinetic Boltzmann equation for the dynamics of large scale in turbulent 
flows? In principle, such an equation should be derived ab-initio through a 
renormalization-group procedure, starting from the true Boltzmann equa- 
tion for molecules. Preliminary attempts in this ambitious direction look 
very encouraging, and yet not conclusive |2!- While waiting for a rigor- 
ous derivation, a practical strategy is to resort to discrete versions of the 
Boltzmann equation, now known as Lattice Boltzmann equation (LBE), 
and endow it with a self-consistent relaxation operator [S|. The crucial 
asset of LBE is that it provides a minimal form of kinetic theory compati- 
ble with the physics of turbulent flows, and such that it can be simulated 
very efficiently on present-day computers. This approach, sometimes called 
LBE-T, has been recently shown to provide leading-edge numerical results 
for turbulent flows in real-life geometries [Sj. Despite the impressive re- 
sults, these simulations leave many theoretical ends open, and 'whether a 
theory can be developed remains to be seen' jH]. Leaving aside the im- 
portant issue of numerical efficiency, from a purely theoretical perspective, 
a relevant question is whether the added value of the kinetic approach to 
fluid turbulence can be linked to the the dynamics of non-hydrodynamic 
fields, high order kinetic moments of the Boltzmann distribution, some- 
times called ghost fields ^1 El H^l EI ■ Ghost fields represent the hidden 
kinetic information which, although necessary to guarantee the correct lo- 
cal symmetries, does not normally emerge to the macroscopic scale. In this 



paper we aim at investigating the role of the dynamics of ghost fields in the 
kinetic approach to fluid turbulence. In particular, we focus our attention 
on the effects of ghost fields on the small-scale statistics of high Reynolds 
number flows. Unfortunately, due to limitations in computational power, 
it is still impossible nowadays to perform numerical simulation in the fully 
developed turbulent regime. To reach high Reynolds numbers flows, with 
a clear separation between the energy injection and the energy dissipation 
scales, one needs to resort to some model of turbulence. A popular class 
of deterministic dynamical models, widely used in recent years, is given by 
"shell-models" ^J . Shell models represent the only example where flows at 
high Reynolds numbers, with realistic small-scales statistics, can be stud- 
ied. They are a good test-bed for new theories and numerical schemes 
aimed at improving the understanding of small-scales turbulent behaviour. 
The paper is organized as follows. First we propose a Lattice Boltzmann 
scheme for shell-models. Second, we study analitycally and numerically 
its hydrodynamical limit with particular emphasis on the importance of 
ghosts fields for the small scales dynamics. Finally, we present numeri- 
cal and analytical results on the "multi-relaxation" regime, i.e. the case 
when hydrodynamic fields and ghost fields have two different relaxation 
properties to the local equilibrium. 

2 Shell models 

Shell models are the simplest dynamical systems featuring a realistic picture 
of energy transfer from large to small scales |15[ [T^ . The main advantage 
of shell models is that they can be analyzed in great depth via highly 
accurate numerical simulations, which permit to resolve up to six decades 
in momentum space. They are non-linear deterministic models of Navier- 
Stokes equations, built such as to preserve energy, helicity and phase-space 
volume in the inviscid and unforced limit. In this work we shall focus on 
the following shell model [T()] : 

dtUn = iknQn " J^fc^C/« + F^ (1) 

where the non linear term is given by: 

Qn = K_^Un-l + ^U*_,Un+l + ^-^-^Un+lUn+2. (2) 

In the above C/„ is a complex variable representing the fluctuating velocity 
field at wavenumber kn — 2"A;o, and & is a free parameter fixing the phys- 
ical dimension of the second inviscid invariant (here we fixed b — —0.5). 
In the presence of a large-scale forcing i^„ = F5n,o, this model exhibits 



excellent scaling laws, from many aspects indistinguishable from those of 
real turbulence. For example, the p-th order structure functions 

5pH = (|c/„n^fc-';(p) (3) 

are characterized by a set of scaling exponents, C{p), very close to those 
measured on turbulent flows. Many other aspects, connected to the velocity 
probability density functions, energy dissipation statistics and multi-time 
multi-scale correlation functions are also in good agreement with what mea- 
sured in experimental and numerical studies of Navier-Stokes equations. 
For these reasons shell models have represented a unique occasion to inves- 
tigate small-scale turbulent statistics without the difficulties of the original 
Navier-Stokes equations. Despite their apparent simplicity, a full system- 
atic analytical control of the small-scale statistical behaviour is still lacking. 
Recently, a series of promising closure attempts based on stochastic closures 
have been proposed {T71 I18L ITOj . Yet, they cannot be considered conclu- 
sive. The importance of Lattice Boltzmann schemes for the shell-model 
||2J) is therefore twofold. First, they may shed some lights on the complex 
multi-time dynamics of the hydrodynamical limit, second they may be use- 
ful to control, and optimize, convergence to the hydrodynamic limit, which 
may be useful also for LB schemes of Navier-Stokes equations. 

3 LBE shell model 

In this section we shall develop a discrete kinetic model whose hydrody- 
namic limit is precisely the shell model eq.Q. To this purpose we introduce 
a 5-speed lattice Boltzmann scheme in the wave number space, kn , obeying 
the following dynamic equations: 

dtfjikn) + tkrjjikn) - --(/,(fc„) - /f (fc„)) (4) 

where fj = [/o, /i, /2, /a, A] is the discrete distribution associated to the 
discrete speeds Cj = [0, 1, —1, 2, —2]. The local equilibrium is given by 

/f (A:„) = wj K + c,Un + (c/ - 1)A.] (5) 

where Wj = [6/12,2/12,2/12,1/12,1/12] are normalized weights. The 
equivalent of the macroscopic density and velocity fields are defined as 
follows: 

4 4 

Rn = E /^(^") Un^Yl fji''n)Cj. (6) 

The third macroscopic field, the analogue of the traceless momentum flux 
tensor, must be adjusted in such a way as to reproduce the non-linear term 



Qn in the shell model. After simple algebra, one derives: 

2D„ + i?„ = Q„. (7) 

Let us notice that the shell- model equations we want to mimick are written 
in terms of a complex variable, C/„. Therefore, here the /j(A;„)'s loose the 
nature of probability density functions they usually have in LB schemes in 
real space (complex distribution functions have been already used in the 
past to simulate quantum mechanics [2^1^ ). 

The first consequence is that in order to keep the "macroscopic density" 
for shell n, i?„, constant in time, one has to modify the streaming term in 
(0J for rest particles as follows: 

5t/o(fc„) = -zA:„(2/3(fc„)-2/4(fc„) + /i(fc„)-/2(fc„))-i(/o(fc„)-/o''(fc„)) 

T 

(8) 
By a linear transformation, we move to the momenta representation for the 
stress tensor, S'„ = J2j fj{f^n)cj, and two ghost fields. An — fi{kn) — fzikn) 
and Bn = fsikn) + fi{kn)- The resulting equations are: 

dtRn = (9) 

dtUn = iknSn (10) 

dtSn = iknUn + GiknAn {Sn - Qn) (H) 

T 

dtAn^2tknBn--{An~l-Un) (12) 

T 3 

dtBn = 2iknAn - -{Bn - ^Qn) (13) 

T 4 

The set of equations H9I13|I is our kinetic shell model. From the first three 
equations, we obtain 

dtUn + TdttUn = iknQn " rklUn " Grfc^ A„. (14) 

This is the master equation of our treatment. First of all, we inspect its 
hydrodynamic limit. To this purpose, we notice that in the limit r — )■ 
the ghost field An collapses onto its attractor, the velocity field: 

limA„ = % (15) 

SO that 1)14(1 delivers: 

dtUn - iknQn - irklUn - rduUn- (16) 



It is therefore seen that, upon neglecting the term rdttUn, which is indeed 
higher order in the Knudsen number, the correct shell model is reproduced 
in the limit Kn -^ 0, with the ghost field contributing a factor 2t to the 
flow viscosity. A perturbative expansion in r of all fields appearing in 
p4|l . reveals that the effect of ghost fields at second order in Kn, yields a 
non-conservative contribution of the form 

T^klQn. (17) 

The finite-Knudsen regime is therefore characterized by the interplay of 
three distinct terms, namely: 

knQn {inertial term) (18) 

rk^JJn {dissipative term) (19) 

T k^Qn {ghost contribution). (20) 

Dimensional matching of these competing terms delivers the relevant crossover 
scales in Fourier space: 

• Dissipative scale (Dissipation=Inertia): 

^ (21) 



r3/4 



• Ghost scale (Ghost term== Inertia): 

fc, ~ i. (22) 

k 

These relations show that kg ~ -^jj >> kd for r << 1 which means that 
the ghost fields cannot play any role on the dissipation properties of the 
system since they do not reach up to the dissipative scale. In order to elicit a 
non-trivial role for the ghost fields, we need to realize the condition kg < k^- 
This necessarily leads to a generalized LBE in which ghost fields relax on 
their own timescale, longer than the hydrodynamic one. The minimal such 
choice is to define two relaxation times: r^ and Tg, for hydrodynamic and 
ghost fields respectively. Since we aim at a fully turbulent regime, we shall 
consider for r^ the smallest possible values compatible with grid resolution. 
The relaxation time for the ghost field will then be changed in order to 
investigate its effects on the dissipation properties of the system. 



4 Multi-relaxation shell BGK model 

The simplest multi-relaxation kinetic shell model takes the following two- 
time form: 

dtRn = (23) 

dtUn = iknSn (24) 

dtSn = iknUn + GiknAn {Sn - Qn) (25) 

dtAn^2lknBn {An ^ -U,,) (26) 

Tg 3 

dtBn^2iknAn (S„--Q„). (27) 

This set of equations is easily reproduced by going back to earliest LB 
formulations, in which coUisional effects were taken into account through a 
scattering matrix Mji describing the interaction between the j'-th and i-th 
populations: 

dtf,{kn) + iknfAkn) = M^,{f,{kn) - f!\kn)). (28) 

Following the top-down procedure introduced in j26| . we can construct a 
scattering matrix with eigenvalues A = {0,0, —1/t^, — I/t^, — 1/rg} and a 
corresponding set of kinetic eigenvectors, V , k — 0,4, associated with 
the set of fields Rn,Un,Sn,An,Bn, respectively: 

/,(fc„) = RnVp + UnVp + SnVp + AnVl'^ + BnV^. 

This corresponds to a partition of the five-dimensional kinetic space into a 
hierarchy of two conserved quantities {Rn, Un), one quasi-conserved (trans- 
port) quantity {Sn) and two ghost fields (A„,B„). Using (|^ and (P5|l we 
obtain: 

dtUn = iknQn " T^k^Un " 6r,fc2 A„ - T^dttUn- (29) 

where the dependence on Tg is implicitly hidden within the fields An and 
Un- It is worth pointing out that at this stage we are still dealing with exact 
equations. In order to get a first guess on the dynamics of ghost fields in 
this case, we make the approximation of imposing steady-state conditions 
on the ghost fields equations H26I27|1 : 

= 2tknBn~—{An-lUn) (30) 

Tg 3 

= 2iknAn--{Bn-^Qn). (31) 



This yields: 

This steady-state approximation must be understood as an estimate for the 
mean value of the fields. We show later, by direct numerical simulations of 
the eqs. (|23l27ll . that the prediction extracted out of (|30|l and (J2U, yields 
the correct qualitative and quantitative statistical behaviours. It is also 
quickly checked that in the limit Tg = t^ << 1 the relation 1)32(1 reduces to 
the expression pSII . as it should. For further analysis it proves convenient 
to explore the behaviour of H32|) in the two regimes of small and large scales 
separately. 

4.1 Large scales 

Setting Tg = 1, large scales are identified by the condition k„ << 1. In this 
regime, the denominator of (|32|l simplifies, P{knTg) -^ 1, and we obtain: 

Upon substituting this in the master equation (|29|) we get: 

dtUn = iknQni^ - iiknTgT^) - Sr^klUn - T^duUn- (34) 

It is therefore apparent that, since t^ << Tg — I and fc„ << 1, the relative 
correction to the convective term, 3TgT^k^, is negligible. As a result, we 
come to the conclusion that the large-scale regime is virtually uncontami- 
nated by the ghost fields. 

4.2 Small scales 

Small scales are identified by the condition A;„ >> 1, again with the position 
Tg = 1. In this regime the denominator P{knTg) is dominated by the kf^ 
term, so that the master equation delivers: 

Upon substituting in (|29|l . we obtain 

r 2r 

dtUn = ikn{l -)Qn ~ T^klUn f C/„ - TdtfUn- (36) 

^9 ^g 

^From this expression we notice that the renormalized convective term is 
still conservative, with a renormalized factor 1 — T^JTg ^ 1. The viscous 



term becomes —Ti^kfJJm corresponding to a viscosity v = Tj^, i.e. three 
times lower than in the previous case. Apart from the usual second order 
time derivative of the velocity field, the remaining piece, — 2r^/Tg[/„, is a 
sub-leading damping term, due to the combined effect of high wavenumbers 
and long ghost relaxation time. In other words, ghosts can act at scales 
larger than the dissipative scale, but their amplitude is suppressed by a 
factor TvJTg and consequently they gently disappear from the scene, leaving 
the system with the bare ghost-free viscosity t^,. Summarizing, the present 
analysis leads to the following predictions: 

• Hydrodynamic scenario: r^ — r^ << 1, ghost fields are enslaved 
to their local equilibrium values. They contribute an extra term 
2rj/ to the fluid viscosity and do not affect the convective terms to 
any appreciable extent. At small but finite Knudsen numbers they 
are confined to sub-dissipative scales only and cannot produce any 
further appreciable effect. As a result, the correct hydrodynamic 
limit is recovered, with an enhanced viscosity v = 3t^. 

• Non-hydrodynamic scenario: t^ — > 0,Tg — 1, ghost fields are no 
longer enslaved to the fiuid velocity. They receive contributions from 
the velocity field and the non-linear term Q„ at all scales through the 
propagator P{knTg). As a result, they exhibit high-frequency, small- 
amplitude fluctuations, which do not affect the large scale behavior 
of the system because they are suppressed by a t^/tq factor. The 
correct hydrodynamic limit is still recovered, with a bare viscosity 
three times smaller than in the the previous case, vq — t^. 



5 Numerical results 

The theoretical scenario discussed in the previous section has been tested 
against numerical simulations of the kinetic shell model. As a first test, 
we have simulated the kinetic shell model in the hydrodynamic regime, 
namely t^ = Tg = 5 x 10^*, with fc„ = 2"~^^, n — 1, 25. For the sake of a 
quantitative comparison, the same simulations have been repeated with the 
original shell model ^ (with a viscosity ly = 3t^). In Figure^we show the 
energy spectra for the two cases. An excellent agreement between the LB 
and the shell model simulations is observed across the whole range of scales, 
except for scales well inside the dissipative range. In the inset of figure 1^ 
we also present a check of the enslaving relation 1)15(1 by plotting the ratio 
between the ghost field and the velocity field, 3i?e(A„)/i?e([/„) for two 
typical wavenumbers, at large scales, n = 6 and at small scales close to the 
dissipative cut-off, n ~ 18. Notice how the relation (|15|l is perfectly verified 
at large scales, while some, small, deviation from slaving is observed at the 



end of the inertial range. This confirms our theoretical analysis, namely 
that ghost fields are completely enslaved to their equilibrium values and 
do not affect the hydrodynamic behaviour of the turbulent system. Only 
strongly dissipative physics is affected by the ghost dynamics. The global 
stability is not changed. 

More interesting, is to explore numerically the non-hydrodynamic regime, 
Tg ^ 0(1). We have performed a set of simulations with, Tj, = 5 x 10~^, 
Tg ^ 0{1). The corresponding spectra for the velocity field are shown in 
Figure 121 where the case t^, = Tg is also reported to highlight the effect of 
the reduced viscosity from i/ = St^, to i/q = t^. As predicted by our theo- 
retical analysis, the LB model in the non-hydrodynamic regime reproduces 
turbulent shell dynamics with the correct hydrodynamic viscosity vq ~ t^, 
while the ghost field tends to decrease in intensity by going to smaller and 
smaller scales. Now ghosts are decoupled from the velocity fluctuations 
and the enslaving relation H15|l is no longer satisfied. Still, they do not 
have any significant impact on the dynamics of the velocity field because 
their intensity is neghgible (see inset of Figure [2Jl. 

A deeper insight into the role of ghost fields at all scales, inertial and dissi- 
pative, can be gained by inspecting the structure functions Q- In partic- 
ular, in figure 131 we show the fourth and sixth order fiatness F„ = -njijj 

and Fn = ttttts as a function of n for the hydrodynamic and non- 
hydrodynamic regime. The shell-model results are also presented. First 
we observe that both quantities have a clear dependency on the scale, that 
is a sign of intermittency in the velocity statistics. Second, all numerical 
results agree in the two regimes, the only difference being an increased 
inertial range extension when ghosts fields decouple. 



6 Conclusions 

Summarizing, we have presented a detailed analysis of the effects of non- 
hydrodynamic (ghost) fields on the statistical properties of hydrodynamic 
turbulence, within the framework of shell models. As a first result, we 
have shown that if ghost fields relax on the same time-scale as the hy- 
drodynamic one, Tg = Ty (the common scenario in current real-space LB 
simulations), then in the hydrodynamic limit where this scale is sent to 
zero, the ghost fields contribute to the hydrodynamic viscosity with a ratio 
2 : 1 with respect to the weight of the hydrodynamic fields: v = 2Tg + Tv 
Higher order ghost contributions, at finite relaxation times, are segregated 
to sub-dissipative scales, so that no further effects on the hydrodynamic 
behaviour can result. The non-hydrodynamic 'overtake' scenario [kg < kd) 
can be realized by allowing ghost fields a longer lifetime than hydrody- 
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namic modes. It is then found that ghost fields do not spoil the inertial 
physics (large scales) up to a neghgible correction r^/rg. The dissipative 
properties of the turbulent system are also reproduced, but now without 
the ghost contribution, leading to an hydrodynamic bare viscosity smaller 
than in the previous case: vq = t^. 

Our analysis supports the counter-intuitive notion that letting ghosts 
'alive' on long time scales leads to a reduction of the viscosity, as compared 
to the case in which ghosts are frozen to their hydrodynamic equilibria. 

Once extrapolated to the framework of real-space lattice Boltzmann 
models, our findings provide a motivation towards multirelaxation models, 
as opposed to the currently popular single-time relaxation BGK model. 
This come-back of multirelaxation models has been invoked by other groups 
as well 123, based on motivations of improved (linear) stability. 

It is tempting to speculate that real-space LB simulations of two and 
three-dimensional turbulence might profit by moving (back) to a multi- 
relaxation scenario in which the hydrodynamic scale is kept to its minimum 
fixed by grid resolution, while the ghost fields time-scale is made much 
longer. 

On a similar vein, one may speculate that making the hydrodynamic 
and ghost time scales respond self-consistently to turbulence observables, 
but with distinct functional dependencies, may prove beneficial also for 
LBGK-r simulations with turbulence modeling [5]. As a final remark, our 
results provide a clear evidence that the hydrodynamic manifold is very 
robust against large fiuctuations of non-hyrodynamic fields. 
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FIGURE 1 
Comparison between the original shell model and our LB model. We 
plot log2(|C/np) vs n for the shell model with viscosity i/ = 3t^ (+) and the 
LB kinetic model with r^y = Tg ^ 5 x 10"'* (x). Both models have the same 
forcing acting on the first two shells. Inset: check of the ghost-velocity 
slaving in the hydrodynamic regime. We plot _Re(3A„(i))/i?e(C/„(t)) vs 
t for shell index n = 6 (straight line) and n = 18 (dotted line) with 
T^ = Tg = 5 X 10^''. Notice the small deviation observed at the small- 
est scale. 

FIGURE 2 
Velocity spectra for hydrodynamic and non-hydrodynamic LB regime. We 
plot log2(|C/„P) vs n for the LB model with the following choice of param- 
eters: T^ — Tg — 5 X 10"'' (-f-) and t^ ~ 5 x 10"'', Tg — 1 (x). Notice the 
increase of the inertial-range extension in the multi-relaxation LB model 
because the ghost field (•) is completely negligible at small scales. Inset: 
check of the ghost-velocity slaving in the non-hydrodynamic regime. We 
plot Re{3An{t)) / Re{Unit)) vs t for shell index n = 6 (straight line) and 
n = 18 (dotted line) with Tg = I. Notice now, at difference from the case of 
fig. n] that the velocity and the ghost fields at small scales are decoupled, 
being the overall intensity of the ghost field negligible. 

FIGURE 3 
Analysis of the flatness F^^^ = {Uf,) / {UlY'"^ for different values of p. We 
plot log2 Fn VS n in our kinetic model for the following choice of parame- 
ters: Ty = Tg — b X 10"'' (x) and Ti, — b x 10""*, Tg = 1 (-I-); to emphasize 
the correct behaviour of the system we also plot the same quantity in the 
case of the original shell model with viscosity vq = Ti, {-k). Inset: the same 
cases but for the sixth order flatness Fn = {U^)/{U^)'^. 
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